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SUMMARY 


The  modified  phase-spece  theory  of  reaction  rate*  has  been  used 
to  predict  the  three-body  recombination  and  dissociation  rate  coefficients 
of  the  diatomic  gas  molecules:  H2,  N2,  02,  F2,  C&2,  Br2 ,  I2,  HF,  HCfc,  CO 
and  NO,  in  the  presence  of  argon  as  a  collision  partner.  The  ability  of 
the  theory  to  quantitatively  predict  and  correlate  both  low  temperature 
recombination  rate  measurements  and  high  temperature  dissociation  rate 
measurements  is  substantial.  The  success  of  the  theory  clearly  illustrates 
the  importance  of  the  weak  attractive  forces  between  the  recombining  atoms 
and  argon  atoms  for  recombination  at  low  temperatures,  the  marked  reduction 
in  the  rates  at  high  temperatures  due  to  nonequilibrium  distributions  in 
the  vibrational  state  populations  of  the  molecules,  and  the  major  contribu¬ 
tions  to  reaction  progress  via  electronically  excited  molecular  states  at 
all  temperatures  for  such  molecules  as  N2  and  CO. 
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1.0  INTRODUCTION 

In  this  paper  we  have  used  the  modified  phase-space  theory 

(1-4) 

of  reaction  rates,  as  developed  by  Keck  and  his  co-workers  for 

three-body  reactions,  to  predict  the  Individual  dissociation  and  re¬ 
combination  rate  coefficients  of  a  number  of  diatomic  gases  diluted 
In  an  argon  heat  bath. 

We  have  adopted  the  phase-space  theory  approach  because  It 
represents  the  most  general  method  of  calculating  the  reaction  rates 
of  atomic  and  molecular  systems  whose  Interaction  can  be  described  by 
the  motion  of  a  representative  point  In  the  classical  phase  space  of 
the  reacting  system.  It  can  be  shown  that  the  more  conventional  classic¬ 
al  reaction  rate  theories  such  as:  unlmolecular  decay  theory, ^  absolute 
reaction  rate  theory/6^  and  available  energy  theory,  ^  are  all  special 
cases  of  the  phase-space  theory.  For  a  general  discussion  of  these  points 
we  refer  to  Keck/1^ 

Provided  that  the  interaction  potential  of  the  three-body  system 
can  be  defined,  the  modified  phase-space  theory  which  we  have  used  here 
enables  calculations  to  be  made  of  dissociation  and  recombination  rate 
coefficients  which  do  not  contain  arbltary  constants  and  undetermined 
"sterlc  factors."  Thus,  when  subsequent  comparisons  are  made  with 
reliable  experimental  measurements,  It  Is  usually  possible  to  ascribe 
any  significant  disagreement  between  the  theory  and  experiment  to  Inadequa¬ 
cies  In  the  assumed  form  of  the  basic  three-body  interaction  potential. 

Such  comparisons  thus  provide  a  method  of  making  quantitative  improvements 
to  the  Interaction  potentials. 


As  we  shall  demonstrate  by  our  comparisons,  we  have  been  able 
to  obtain  substantial  quantitative  agreement  between  theory  and  experi¬ 
ment  and  have  thus  correlated  low  temperature  three-body  recombination 
rate  measurements  (obtained  for  example,  using  the  "discharge-flow-tube" 
technique)  and  high  temperature  "shock-tube"  dissociation  rate  measure¬ 
ments  for  those  cases  where  both  types  of  data  are  available.  The  overall 
success  of  the  theory  clearly  illustrates  the  importance  of:  1)  the 
attractive  Van  der  Waal's  forces  between  the  recombining  atoms  and  the 
argon  atoms  for  recombination  at  low  temperatures  (such  forces  are,  of 
course,  invoked  in  the  "relaxed  complex"  mechanism  of  recombination ^  ^), 

2)  the  departure  from  an  equilibrium  distribution  of  the  vibrational  state 
(3  12) 

populations  '  which,  for  either  dissociation  or  recombination,  occurs 

at  an  energy  level  on  the  order  of  kT  below  the  dissociation  limit,  3)  the 

possible  major  contributions  to  the  overall  reaction  rate  due  to  reaction 

(4) 

progress  via  electronically  excited  molecular  states. 

The  source  experimental  measurements  which  we  have  used  in  our 
comparisons  Were  mainly  suggested  by  recent  review  articles.  ^  ^  The 
particular  cases  which  wa  have  considered  are  the  homonuclear  diatomic 
gases:  ,  N^,  0 2>  F^,  Cf^,  Br2,  I^,  and  the  heteronuclear  gases:  HF,  HCl, 

CO  and  NO. 

In  this  paper  we  shall  not  include  the  detailed  development  of  the 

(1-4) 

modified  phase-space  theory  since  this  can  be  found  elsewhere.  A 

discussion  of  the  theory  and  a  complete  summary  of  the  working  formulae 

(4) 

has  been  given  by  Shui,  Appleton  and  Keck  for  the  particular  case  of 
the  dissociation  and  recombination  of  nitrogen  in  an  argon  heat  bath. 
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We  included  in  that  paper  a  general  account  of  the  manner  in  which 
reaction  progress  via  electronically  excited  states  should  be  treated, 
together  with  a  corrected  version  of  the  way  in  which  departures  from 
an  equilibrium  distribution  in  the  vibrational  state  populations  are 
allowed  for.  Thus,  in  Section  2.0  we  shall  give  a  brief  summary  of  the 
working  formulae  and  their  physical  interpretation,  and  devote  the  re¬ 
mainder  of  the  paper  to  discussion  of  the  comparisons  between  theory  and 
experiment  (Section  3.0),  and  our  final  conclusions  (Section  4.0). 
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2.0  THEORETICAL  CONSIDERATIONS 

The  type  of  reactions  with  which  we  are  concerned  are  commonly 
written: 


XY  +  M  T — —  X  +  Y  +  M 
k 


(2.1) 


and  are  observed  to  proceed  in  accordance  with  the  overall  rate  equation: 


d[XY]/dt  -  -  kd[XY](M]  +  kr[X][Y][M] 

-  -  d[X]/dt  -  -d[Y]/dt,  (2.2) 


where  kd  and  kf  are,  respectively,  the  overall  dissociation  and  recombina¬ 
tion  rate  coefficients.  For  hotnonuclear  diatomic  molecules:  X2,  Eq.(2.2) 
becomes: 

d(X2J/dt  -  -  kd[X2][M]  +  kr(X]2[M] 


-  -  d[X]/2dt  . 


(2.3) 


It  is  implied  by  the  above  rate  equations  that  we  restrict  our  attention 

to  the  calculation  of  reaction  rates  at  times  which  are  long  by  comparison 

with  any  "induction  periods,"  i.e.  the  time  Intervals  required  to  achieve 

close  approximations  to  steady  state  distributions  in  the  internal  energy 

states  of  the  molecules  and  atoms  (see,  for  example,  reference  (17)  for  a 

discussion  of  the  induction  period  due  to  vibrational  relaxation).  Therefore, 

both  k.  and  k  are  to  be  regarded  as  functions  of  temperature  alone,  and 
d  r 
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their  ratio  equal  to  an  equilibrium  constant.  If,  as  we  shall  assume 
(4) 

here,  each  allowable  electronic  state  of  a  particular  species  is 
present  in  the  reacting  system  in  its  local  equilibrium  proportion  re¬ 
lative  to  the  total  concentration  of  that  species,  then 

k./k  -  EQ(p)  £Q(q)/£0(m) 
dr  p  q  m 

-  K  (T)  (2.4) 

e 

where  Q(p)  is  the  total  particle  partition  function  for  atom  Xp  in 
electronic  state  p,  and  Q(q)  and  Q(m)  are  similarly  defined  for  atom 
Y  and  molecule  XY  ,  respectively.  The  overall  rate  coefficients  are 

C]  01 

thus  written  in  the  forms: 

kd  "  m  kd*mlp,q)  Q(»)  (2.5) 

and  k  -  I  k  (m|p,q)  Q(p)Q(q)/E  Q(p)E  0(q)  (2.6) 

r  tn  r  p  cj 

where  kd(m|p,q)  is  the  rate  coefficient  for  dissociation  of  the  molecule 
in  electronic  state  m,  to  atoms  in  electronic  states  p  and  q,  and  kr(m|p,q) 
is  the  corresponding  recombination  rate  coefficient.  We  note  that  as  a 
consequence  of  the  Bom-Oppenhelmer  separation,  the  specification  of  a 
particular  molecular  state  automatically  identifies  the  atomic  states  to 
which  the  molecule  dissociates. 

For  most  diatomic  molecules  and  their  atoms  it  is  usually  true 
that  the  ground  state  partition  functions,  e.g.  Q(p-O),  are  much  greater 
than  the  partition  functions  of  the  first  and  higher  excited  states  for 
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temperature  ranges  of  normal  interest,  and  thus,  the  overall  rate 
coefficients  may  be  approximately  rewritten  in  the  forms: 


kd  " 

Z  kd(mj0,0)Q(m)/Q(m-0) 

(2.7) 

and 

k  - 

r 

Z  k  (m|0,0), 
in  r 

(2.8) 

so  that. 

Vs. 

-  Q(p-0)Q(q-0)/Q(m-0) 

(2.9) 

is  the  equilibrium  constant  based  on  ground  state  concentrations  only. 

One  important  exception  for  which  the  above  approximate  equations  are 

not  valid  over  the  whole  temperature  range  of  our  interest  is  the  case 
3 

of  oxygen.  The  P  "ground"  state  oxygen  atom  actually  consists  of  three 

•i  n  n 

states  ?2,  and  PQ,  with  slightly  different  energy  levels:  0,  229°K, 

and  326°K,  and  with  corresponding  degeneracies  of  5,  3,  and  1,  respective¬ 
ly.  Thus,  at  temperatures  below  about  1000°K,  the  exact  equations  (2.5) 
and  (2.6)  have  been  used  here. 

The  modified  phase-space  theory  enables  us  to  calculate  the 
individual  recombination  (or  dissociation)  rate  coefficients:  kr(m|p,q), 
in  terms  of  what  Keck^  has  called  the  "barrier  rate"  coefficient: 
ky  (m|p,q).  This  latter  quantity,  which  represents  a  rigorous  upper  bound 
to  the  actual  rate  coefficient,  requires  modification  by  a  statistical 
correction  factor:  (N/NQ),  and  nonequilibrium  correction  factors:  (k/k^), 
to  obtain  kr(m|p,q). 

The  barrier  rate  is  given  in  terms  of  the  total  one-way  flow  of 
points,  representative  of  the  three-body  complex  (X-Y-M) ,  in  the  classical 


phase  space  across  a  suitably  defined  surface  which  separates  the 
initial  states  of  the  system  -  the  reactants,  from  their  final  states  - 
the  products.  The  phase  space  surface  which  most  logically  separates 
molecules  from  atoms  is  the  "barrier  surface"  which  is  defined  by  the 
constraints 


exy  -  bxy  "  0  ;  E  >  0  (2*11) 

where  E  is  the  total  Internal  energy  of  the  three-body  system  (X-Y-M), 

Exy  Is  the  Internal  energy  of  the  recombining  atoms  (X-Y) ,  and  is 

the  height  of  the  rotational  energy  barrier.  For  recombination,  the 

product  state  is  defined  by  (EXY  -  BXY)<  an<*  the  reactants  bY 

(1  2) 

(Exy  -  Bxy)>  0.  By  following  this  prescription,  Keck  *  has  derived 
the  following  expression  for  the  barrier  rate  coefficient: 

ky(m|p.q)  ■  fa2*2^z2  "  ^MSkT/^yxY*1^2!1  "  exP(“  Bm/kT)  ]  (2.12) 

where  f  -  g^  / gx  g  ,  is  the  electronic  degeneracy  factor,  is  the 
®  p  q 

reduced  mass  of  the  molecule  XY,  and  B  is  the  maximum  height  of  the 

in 

rotational  barrier. 

For  the  purpose  of  providing  a  simple  physical  explanation  for 

2 

the  form  of  Eq.  (2.12),  we  may  point  out  that  the  factor:  Airz^a^-  *j_)  * 

Is  a  molecular  volume  proportional  to  the  number  of  atom  pairs  (X-Y), 

2  1/2 

close  enough  to  recombine;  whereas,  the  factor:  'Ta  (8kT/tryXY)  ,  Is  a 
rate  constant  proportional  to  the  frequency  at  which  the  (X-Y)  pairs  are 
stabilized  under  the  influence  of  the  third  body  M.  The  reason  for  the 


-  fl¬ 


ap  pe  trance  of  Che  reduced  mass  U^,  for  Che  recombining  atoms  In  the 

velocity  term  rather  than  the  reduced  mass  l^.yj'pr  c°lli8ion8  of  (X-Y) 

with  M,  is  that  in  the  "barrier  rate"  it  is  the  rate  of  momentum  transfer 

from  M  to  (X-Y),  rather  than  the  collision  rate  of  M  with  (X-Y),  which 

controls  the  recombination.  The  additional  factor:  [1  -  exp(-  B  /kT)], 

in 

simply  eliminates  those  atom  pairs  which  cannot  form  bound  molecules 
because  of  their  excessive  orbital  angular  momentum. 

The  characteristic  lengths  a,  z^ ,  and  are  all  temperature 
dependent  and  may  be  determined  in  terms  of  the  characteristic  potential 
parameters  which  are  used  to  describe  the  interaction  potential  of  the 
three-body  system  (X-Y-M) .  We  have  assumed  that  the  interaction  potent¬ 
ial  is  given  as  the  sum  of  two  potentials: 

V0  -  VXY  +  V1M  (2’13> 

where  is  the  potential  of  the  recombining  atoms  (X-Y),  and  V^M  is 
the  potential  between  one  of  the  atoms  and  the  third  body  M;  ViM  =  V^, 
when  the  internuclear  distance  r^  is  smaller  than  r^,  and  ViM  =  V^, 
when  r^  >  r^.  Keck  has  termed  this  the  "dumb-bell"  model  of  the  inter¬ 
action  potential.  In  our  calculations  we  have  used  the  Morse  potential 
function: 


V(r)  -  U{[1  -  exp(-  B(r  -  rft))]2  -  1}  (2.14) 


to  represent  both  V^y  and  V^,  with  appropriate  suffices  added  to  the 
potential  parameters  U,  B,  and  r^  to  identify  the  particular  two-body 


irniamwM 
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potentlal  under  consideration.  Graphical  solutions  for  a,  z^,  and  z^ 

as  functions  of  the  appropriate  potential  parameters  have  been  given 
(4) 

elsewhere.  However,  we  should  point  out  that  the  phase-space  theory, 

In  the  form  which  we  have  used,  tacitly  assumes  that  all  species,  includ¬ 
ing  the  complexes:  X.M  and  Y.M,  are  present  in  their  local  equilibrium 
proportions,  and  as  a  consequence,  the  major  temperature  dependence  of 
the  effective  collision  cross-section  (ffa  ),  is  determined  by  a  factor: 
exp(U1M/kT). 

Unfortunately,  the  barrier  surface  which  is  chosen  to  separate 
reactants  from  products  in  the  phase  space  is  not  unique  in  the  sense 
that  representative  point  trajectories  may  cross  and  re-cross  it  several 

g 

times.  Since  is  calculated  on  the  basis  of  the  total  one-way  crossing 
rate,  it  does  represent  an  upper  bound  to  the  actual  rate. 

(18) 

Keck  has  investigated  this  effect  using  Monte  Carlo  methods 
which  sample  phase-space  trajectories  which  cross  the  "barrier"  surface, 
and  by  integration  of  the  classical  equations  of  motion  in  both  time- 
wise  directions,  has  obtained  a  statistical  determination  of  the  fraction 
(N/Nq),  in  which  N  is  the  number  trajectories  which  result  in  a  complete 
one-way  reaction,  and  Nq  is  the  total  number  of  trajectories  sampled. 

Thus,  we  regard  (N/Nq)  as  a  statistical  correction  factor  which  multiplies 

g 

the  value  of  k^  given  by  Eq.  (2.12).  Numerical  values  for  (N/Nq)  can  be 
found  in  reference  (18)  and  semi -empirical  formulae  in  references  (1)  and 
(4). 

The  nonequilibrium  correction  factor  which  we  have  used  here  is 

(3) 

based  on  the  results  of  the  investigation  by  Keck  and  Carrier  of  the 
coupled  vibration-dissociation-recombination  process  for  a  dilute  mixture 
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of  diatomic  molecules  In  a  heat  bath  of  inert  collision  partners. 

After  the  initial  vibrational  relaxation  transient,  it  was  found  that 
for  dissociation,  the  steady-state  vibrational  level  population  is  very 
nearly  Boltzmann  except  near  the  dissociation  limit  (within  a  few  kT  of 
the  dissociation  limit)  where  the  levels  are  underpopulated;  for  recombi¬ 
nation,  the  steady-state  distribution  is  Boltzmann  near  the  dissociation 
limit,  but  the  lower  levels  are  underpopulated.  Because  of  this,  the 
steady-state  rate  coefficients:  and  k^,  are  smaller  than  the  equi¬ 

librium  rate  coefficients:  kde  and  kfe,  although  their  ratio  is  very 
nearly  equal  to  the  equilibrium  constant;  i.e., 

kd/kr  "  kde/kre  "  Ke(T)  <2.15) 

By  assuming  a  classical  Morse  oscillator  to  represent  the  molecule  and 
an  exponential  form  to  represent  the  repulsive  Interaction  between  the 
collision  partner  and  the  molecule,  Keck  and  Carrier  solved  the  appropri¬ 
ate  master  equation  to  give  a  closed  form  expression  for  the  nonequilib¬ 
rium  correction  factors:  (k/k^). 

In  reference  (A)  it  was  pointed  out  that  when  an  attractive  minimum 
exists  in  the  two-body  potentials  V^,  there  are  four  contributions  to  the 
barrier  rate  coefficient,  viz. 

kj(m|p,q)  -  k®x+  +  k®x_  +  k®y+  +  k®y_  (2.16) 

corresponding  to  four  distinctly  different  configurations  of  the  three- 
body  complex  for  which  the  momentum  transfer  rate  between  (X-Y)  and  M 
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are  maximum.  It  is  therefore  necessary  to  evaluate  a  nonequilibrium 
correction  factor,  i.e.  (k/ke)^+,  for  each  of  these  configurations, 
thereby  obtaining  an  expression  for  the  actual  rate  coefficient  in 
the  form: 


kr(»|p,q>  -  [*Jx+(k/ke)X+  +  krX-(k/ke}X- 

(2.17) 

+  k®Y+(k/ke)y+  +  kJY_(k/ke)Y_]  (N/Nq) 
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3.0  COMPARISONS  BETWEEN  THEORY  AND  EXPERIMENT 

The  potential  parameters  r^M,  for  the  Morse  potential  which 

describe  the  interactions  between  the  separated  atoms  X  and  Y,  and 

the  collision  partner  M,  were  first  estimated  from  empirical  rules, 

(19) 

suggested  by  Bernstein  and  Muckerman,  which  related  r^  to  the  radii 

of  the  principal  maxima  in  the  radial  distribution  functions  of  the 

ground  state  atoms.  The  potential  energy  well  depths:  U^M,  were 

(19) 

first  estimated  by  applying  the  Badger-Johnston  relation  in  the 
manner  suggested  by  Bernstein  and  Muckerman.  The  remaining  Morse  potent¬ 
ial  parameters,  were  obtained  by  setting  2(5  ■  l/LiM,  where  is 

the  range  of  an  effective  exponential  potential,  given  by  the  Mason- 

(21) 

Vanderslice  extrapolation  formula.  These  estimates  of  r  U.u  and 

•IN  IN 

8,u  are  summarized  in  Table  I. 

IN 

The  spectroscopic  data  for  the  molecules  which  are  required  in 
the  calculations  are  summarized  in  Table  II.  These  data  were  mostly 
taken  from  Herzberg,^^  Gilmore,  Krupenie,^^  and  the  JANAF  Tables. 

The  source  references  for  the  experimental  data  have,  to  a  large 
extent,  been  taken  from  the  review  articles  cited  previously,  although, 
we  have  attempted  to  Include  in  our  comparisons  some  more  recent  investiga¬ 
tions.  Because  of  space  limitations,  we  clearly  cannot  present  a  critical 
review  of  the  many  experimental  investigations. 

For  those  molecules  where  both  low  temperature  recombination  rate 
data  and  high  temperature  dissociation  rate  data  are  available,  we  have 
made  the  comparisons  on  the  basis  of  the  Inferred  recombination  rate 
coefficients.  For  those  molecules  where  only  dissociation  rate  data  are 
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avallable,  we  have  made  the  comparisons  using  the  usual  Arrhenius  plots. 

Only  for  the  case  of  oxygen  have  both  methods  of  comparison  been  used. 

We  again  stress  that  the  theoretical  predictions  were  made  on  the  assump¬ 
tion  that  argon  was  the  collision  partner  which,  in  the  majority  of  the 
experimental  investigations  chosen  for  comparison,  was  the  case. 

3.1  Hydrogen 

1+12  2 

Figure  1  shows  our  prediction  of  k^CX  E^|  S,  S),  given  by  the 

lower  dashed  line,  as  compared  with  the  direct  recombination  rate  measure- 
(26  27) 

meats v  ’  '  and  the  recombination  rate  coefficients  Inferred  from  the 

( 28-31) 

high  temperature  dissociation  rate  measurements  for  hydrogen.  ' 

The  disagreement  between  our  prediction  and  the  highest  high  temperature 

measurements  by  Hurle,  Jones,  and  Rosenfeldv  '  is  about  a  factor  of  two. 

However,  our  agreement  with  the  experiment  results  of  Sutton,  and 
(31) 

Myerson  and  Watt  and  the  low  temperature  recombination  measurements 

(26  27) 

of  Larkin  and  Thrush v  ‘  is  excellent.  Thus,  in  view  of  the  possible 
experimental  scatter,  as  Indicated  by  the  cross-hatched  area  which  corres¬ 
ponds  to  the  error  estimate  given  in  reference  (28) ,  we  conclude  that  our 
prediction  is  good. 

The  upper  dashed  line  in  Fig.  1  shows  our  prediction  of  the  equi¬ 
librium  rate  coefficient:  k®(X  *E*|^S,^S) (N/Nq) .  By  comparing  this  with 
the  actual  rate  we  see  that  the  correction  due  to  the  nonequilibrium  distri¬ 
bution  in  the  vibrational  state  populations  is  significant  over  the  entire 
temperature  range  of  the  experiments,  but  becomes  most  pronounced  at  high 
temperatures.  The  magnitude  of  the  nonequilibrium  correction  was  found 
to  be  similar  for  all  of  the  homonuclear  molecules  considered,  although 
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we  have  not  made  the  comparisons  In  the  following  figures  In  order  to 
avoid  confusion. 


3.2  Nitrogen 

We  have  treated  the  case  of  nitrogen  In  some  considerable  detail 
(4) 

In  a  previous  paper.  In  that  paper  we  discussed  the  experimental 
(32) 

evidence  which  Indicated  that  direct  recombination  to  the  first 
excited  state  of  the  molecule  exceeded  that  to  the  ground  state  by  at 
least  50%,  and  also  Included  a  discussion  of  the  kinetic  mechanism  where¬ 
by  the  first  excited  state  Is  populated  and  depopulated  by  collisions. 
Figure  2  Is  taken  directly  from  reference  (4).  The  upper  dashed  curve 
Is  our  prediction  of:  [k^. (X  *Z*|^S,^S)  +  k^(A  ^Z*|^S,^S)],  and  the  lower 

dashed  curve  Is  the  prediction  of  k  (X  1Z+|^S,^S),  alone.  Both  of  these 

r  g 

curves  were  calculated  using  the  empirically  derived  potential  parameters 

given  In  Table  I.  We  have  calculated  the  contributions  to  the  overall 

recombination  rate  coefficient  due  to  recombination  to  other  excited 
5  +  3 

states,  e.g.,  £  ,  A  ,  etc.,  and  have  found  them  to  be  negligible  over 

8  u 

the  entire  temperature  range  considered,  l.e.,  90°K  to  20,000°K. 


Since  the  temperature  dependence  of  our  prediction  of  k^  does 

(32) 

not  match  that  of  the  measurements  due  to  Campbell  and  Thrush  and 
(33) 

Clyne  and  Stedman  particularly  well,  and  In  view  of  our  previous 

comment  that  the  phase  space  theory  In  the  form  used  shows  that  the 

temperature  dependence  of  kr  at  low  temperatures  Is  primarily  determined 

by  the  factor:  exp(U„.  /kT) ,  we  Increased  the  value  of  U...  (see  Table  I) 

and  thereby  obtained  the  predictions  Illustrated  by  the  full  curves  In 

Fig.  2.  By  making  changes  In  both  of  the  other  two  parameters  l.e. 

r  ...  and  0...  ,  theoretical  values  of  k  can  be  obtained  which  correlate 
eNAr  NAr  r 
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any  given  set  of  low  and  high  temperature  data.  Our  preference  for 

(A) 

only  changing  U„.  ,  resulted  from  our  observation  that  the  high 

ilAL 

temperature  prediction  of  thus  obtained  agreed  remarkably  well  with 

(36) 

the  measurements  due  to  Appleton,  Steinberg  and  Llquomik,  and  be¬ 
cause  there  is  recent  Independent  evidence  from  molecular  beam  scattering 

(37) 

experiments  due  to  Jordan,  Colgate,  Amdur  and  Mason,  which  supports 

(21) 

the  value  of  given  by  the  Mason  Vandersllce  extrapolation  formula. 

3.3  Oxygen 

Figure  3  shows  our  prediction  of  ky(X  P) ,  lower  dashed 

curve,  and  our  prediction  of  the  sum  of  the  individual  recombination 
rate  coefficients  corresponding  to  all  the  electronic  states  of  oxygen 
listed  in  Table  II,  upper  dashed  curve. 


Improved  agreement  with  the  temperature  dependence  of  the  low 

(32) 

temperature  measurements  due  to  Campbell  and  Thrush  could  be  obtained 
by  increasing  our  estimate  of  UQAr.  However,  this  would  also  entail  a 
reduction  in  our  estimate  of  r egAr  and  a  change  in  3QAr  for  the  theoreti¬ 
cal  value  of  kf  to  still  pass  through  the  main  body  of  the  high  tempera¬ 
ture  data.  We  have  not  pursued  this  curve  fitting  procedure  in  the  case 
of  oxygen  since  we  have  rather  arbitarily  assumed  that  each  combination 

3 

of  the  P.  .  _  states  of  the  oxygen  atoms  Interact  without  restriction 
/  9 1  tu 

along  each  of  the  potential  energy  curves  of  the  molecule  listed  in 
Table  II.  We  have  thus  violated  the  "non-crossing"  rule  of  potential 
energy  surfaces  which  may  result  in  significant  errors  in  our  calculation 


of  kr  at  the  low  temperatures  where  the  energy  spacing  between  the  P 


states  are  comparable  with  kT 
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It  is  Interesting  to  note  that  the  recombination  rates  Inferred 
from  the  dissociation  rate  measurements  due  to  Watt  and  Myerson^^ 
agree  quite  well  with  the  prediction  of  kr(JE~|  r,  P) ,  lower  dashed 
curve.  This  is  indeed  precisely  what  might  be  expected,  since  Watt 
and  Myerson  used  the  technique  of  atomic  line  resonance  absorption  to 
monitor  the  production  of  oxygen  atoms  behind  shock  waves.  Due  to  the 
high  sensitivity  of  the  technique,  they  were  only  able  to  measure  the 
initial  rate  of  dissociation  over  a  period  which,  although  encompassing 
the  ground  state  vibrational  relaxation  transient,  quite  possibly  did 
not  extend  into  the  region  where  the  higher  electronic  states  assumed 
their  local  equilibrium  concentrations. 

Figure  4  shows  an  Arrhenius  plot  of  the  various  dissociation  rate 
(40-43) 

measurements  together  with  our  theoretical  prediction  of  k.. 

a 

(41) 

It  is  apparent  that  the  theory  underestimates  some  of  Wray's  measure- 

(42) 

ments  and  Camac  and  Vaughan's  measurements  by  about  a  factor  of  two 
to  three.  However,  it  should  be  pointed  out  that  within  the  temperature 
range  where  the  measurements  overlap,  there  is  disagreement  between  them 
of  at  least  a  factor  of  two. 

3.4  Fluorine .  Chlorine .  Bromine ,  and  Iodine 

An  Arrhenius  plot  of  the  dissociation  rate  measurements  for  fluorine 

(44) 

due  to  Johnson  and  Britton  is  shown  in  Figure  4  together  with  our  pre- 
^  +  2  2 

diction  of  kj(X  P,  P) .  Although  there  is  a  relatively  large  scatter 

amongst  the  experimental  data,  the  theoretical  prediction  does  pass  through 
the  main  body  of  the  data  and  thus  appears  to  be  satisfactory. 

The  comparisons  between  the  a  priori  theoretical  predictions  of 
the  recombination  rate  coefficients  (i.e.  theoretical  predictions  obtained 
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uslng  potential  parameters  derived  from  empirical  rules  -  dashed  lines 
in  Figs.  6-8)  and  the  experimental  measurements  for  Ci^t  Br2,  I2, 
illustrate  a  definite  trend.  For  the  case  of  chlorine,  Fig.  6,  the 
theoretical  prediction  lies  below  the  experimental  data  over  the  major 
portion  of  the  temperature  range  and  has  a  less  steep  temperature  depend¬ 
ence  than  that  indicated  by  the  combined  measurements.  On  the  other  hand, 

1  +  2  2 

Fig.  8  for  iodine  shows  that  the  a  priori  prediction  of  [kr(X  £  |  P^2*  ^3/2 
3  2  2 

k^(A  niu|  P3/2*  P3/2^  8reater  than  the  measurements  by  about  a  factor 
of  two,  but  exhibits  a  temperature  dependence  which  is  in  reasonable  agree¬ 
ment  with  the  combined  measurements.  The  comparison  for  the  case  of  bromine, 
Fig.  7,  appears  to  be  an  intermediate  situation,  l.e.  the  prediction  of  the 
recombination  rate  is  greater  than  the  hulk  of  the  experimental  data  at 
temperatures  above  300°K,  but  at  lower  temperatures  it  lies  below  an  extra¬ 
polation  of  the  experimental  results. 

Having  recognised  these  trends  and  again  in  view  of  our  previous 

comment  that  the  major  temperature  dependence  of  ky  is  largely  determined 

by  the  factor:  exp(UiM/kT),  at  low  temperatures,  whereas,  it  is  roughly 
2 

proportional  to  re^M»  independently  of  temperature,  we  adjusted  the  values 
of  U1M  and  rfilM  (i  -  Cfc,  Br,  I)  as  shown  in  Table  I,  to  obtain  the  full 
theoretical  curves  illustrated  in  Figs.  6-8.  It  is  apparent  that  this 
fitting  procedure  has  enabled  us  to  correlate  a  large  fraction  of  the 
experimental  data  for  dissociation  and  recombination  of  the  halogens, 
and  as  a  consequence,  we  suggest  that  the  derived  potential  parameters, 
particularly  and  UgrAr>  are  more  realistic  than  those  given  by  the 


empirical  rules. 
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The  above  curve  fitting  procedure  follows  In  much  the  same  spirit 
as  Porter chose  to  view  the  binding  energy  of  the  relaxed  complex: 

I.M,  as  an  adjustable  parameter  for  correlating  iodine  recombination 
measurements,  and  more  recently,  as  Ip  and  Burns similarly  chose 
to  view  the  binding  energy  of  the  relaxed  complex:  Br.M,  as  an  adjust¬ 
able  parameter  for  correlating  their  bromine  recombination  measurements. 
For  the  I.Ar  complex,  Porter  suggested:  U_  .  a  500°K,  and  for  Br.Ar,  Ip 

1  e  AIT 

and  Burns  suggested:  Ugr  ^  -  A50°K,  c.f.  Table  I. 

3.5  Hydrogen  Fluoride  and  Hydrogen  Chloride 

The  Arrhenius  plots  of  the  dissociation  rates  of  hydrogen  fluoride 
and  hydrogen  chloride,  illustrated  in  Figures  9  and  10,  show  that  our 
predictions  of  the  dissociation  rate  coefficients  are  much  too  low  by 
comparison  with  the  experimental  measurements.  Wg  have>  however, 

also  plotted  the  barrier  rate  coefficients:  kj(X  £  I  S,  P) ,  which,  al¬ 
though  they  are  too  large  (as  we  should  expect),  nevertheless,  have  about 
the  right  temperature  dependence.  We  believe  that  the  reason  why  our 
predictions  of  the  rate  coefficients  underestimate  the  measured  rates  is 
that  our  estimate  of  the  statistical  correction  factors:  (N/NQ)  <  0.01, 
and  the  nonequilibrium  correction  factors:  (k/k^)  —  0.01,  are  too  small. 
The  correlation  formulae  which  we  have  to  estimate  the  values 

of  (N/Nq)  and  (k/kg)  were  deduced  using  the  results  of  Monte  Carlo  trajec¬ 
tory  calculations  and  master  equation  solutions  in  which  the  recombining 
molecules  and  the  collision  partners  had  similar  masses.  We  believe  that 
the  mechanics  of  reacting  collision  processes  in  which  one  of  the  recom¬ 
bining  atoms  has  a  very  small  mass  is  sufficiently  different  from  those 
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cases  considered  in  reference  08),  that  the  correlation  formula  Is 
no  longer  valid.  Further  work  on  this  problem  Is  In  progress. 

3.6  Carbon  Monoxide  and  Nitric  Oxide 

Figure  11  shows  an  Arrhenius  plot  of  the  measured  dissociation 

rate  coefficients (66-69)  carbon  monoxide  together  with  our  predictions 

of  k.(X  *Z+| ^P,^P) ,  lower  dashed  curve,  and  Z  k(m|^P,^P),  upper  dashed 
a  g  m 

curve  (see  Table  II  for  identification  of  the  electronic  states  included). 
The  prediction  which  includes  all  of  the  excited  state  contributions  is 
nearly  an  order  of  magnitude  greater  than  k^(X  Z^|  P,  P)  and,  consequent¬ 
ly,  is  in  better  agreement  with  the  measurements.  Small  increases  in 

either  r  ...  or  could  be  made  which  would  allow  the  theoretical  value 
eiM  iM 

3  3 

of  Z  k.(m|  P,  P)  to  correlate  quite  well  with  the  bulk  of  the  experimental 
tn  d 

<5  o 

data.  The  fact  that  Z  k.(m|  P,  P)  does  agree  so  well  with  the  experiment- 

m  a 

/r  q\ 

al  measurements  is  quite  gratifying  since  Fairbainr  ‘'and  Appleton, 
Steinberg,  and  Liquomik^^  (point  values  in  Fig.  11)  did  observe  an 
induction  time  immediately  following  the  shock  wave,  prior  to  the  onset 
of  a  steady-state  dissociation,  which  they  attributed  to  the  time  required 
to  populate  excited  molecular  states. 

Figure  12  shows  our  predictions  of  k^(X  ^l|^S,^P),  lower  dashed 

curve,  and  [k^(X  ^Ilfs,^P)  +  k^(a  ^Jl|^S,^P)],  upper  dashed  curve,  together 

with  the  few  experimental  rate  measurements  for  nitric  oxide  which  are 
(38  70) 

available.  ’  '  The  agreement  is  fairly  good,  although,  it  would 
certainly  be  desirable  to  have  more  extensive  measurements  of  the  rates 
of  this  important  molecule.  Again,  as  in  the  case  of  CO,  a  slight  in¬ 
crease  in  the  value  of  would  raise  the  high  temperature  rate  coeffi¬ 
cients  and  place  the  theoretical  predictions  in  closer  agreement  with  Wray's 
dissociation  rate  measurements  . 
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4.0  CONCLUSIONS 

In  the  foregoing,  we  have  demonstrated  that  modified  phase-space 
trajectory  calculations  of  the  reaction  rate  coefficients  of  diatomic 
molecules  in  an  argon  heat  bath  are  able  to  correlate  both  the  very 
high  and  low  temperature  rate  measurements.  The  exceptions  to  this 
conclusion  appear  to  be  those  diatomic  molecules  which  have  disparate 
atomic  masses  such  as  HF  and  HCl. 

Apart  from  the  recognition  that  reaction  progress  via  excited 
molecular  states  is  likely  to  be  important  in  certain  cases,  e.g.  N^ 
and  CO,  the  quantitative  success  of  the  theory  is  clearly  dependent  on 
prior  knowledge  of  the  interaction  potential  between  the  recombining 
atom  pairs  and  the  collision  partner.  However,  on  the  basis  of  the 
comparisons  presented  here,  we  suggest  that  valid  quantitative  informa¬ 
tion  on  the  Interatomic  potentials  can  be  obtained  from  a  proper  synthesis 
of  the  phase-space  theory  and  reliable  experimental  data  obtained  over  a 
wide  temperature  range.  Such  a  view  is,  of  course,  entirely  in  keeping 
with  existing  methods  which  have  been  widely  used  for  deducing  inter- 
molecular  potential  information  at  thermal  energies  from  measurements 
of  transport  properties. 

Finally,  we  should  point  out  that  for  those  cases  where  the 
collision  partners  are  highly  reactive,  e.g. 

X  +  X  +  X  ?=^X2  +  X, 

the 'Bumb-bell"  model  for  the  three-body  Interaction  potential  and  the 

£ 

approximations  used  in  deriving  the  expression  for  kf  are  no  longer 
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valld.  A  more  realistic  model  for  the  interaction  potential  must  be 
assumed  and  a  more  detailed  analysis  will  be  involved  in  the  calculation. 
Work  in  this  direction  is  being  persued  at  this  laboratory  on  the 
(H  +  H  +  H)  system  which,  hopefully,  will  shed  some  light  on  the  question 
as  to  why,  in  dissociation  and  recombination  processes,  the  parent  atoms 
are  much  more  efficient  third-body  collision  partners  than  those  which 


are  "inert". 
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FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5 


Comparison  between  theoretical  predictions  of  the  recombina¬ 
tion  rate  coefficient  and  experimental  measurements  for  Hg, 

(26)  (27) 

as  a  function  of  T.  0,  Larkin  and  Thrush/  '  and  Larkin; 

(H),  Hurle  et  al.(28)  (J) ,  Jacobs  et  al;(29)  (M;S),  Myerson 

and  Watt/31^  and  Sutton. ^30^  Dashed  curves:  see  text,  Sec.  3.1. 

Comparison  between  theoretical  predictions  of  the  recombination 

rate  coefficient  and  experimental  measure-ments  for  N^,  as  a 

(33) 

function  of  T.  0,  Clyne  and  Stedman;  f,  Campbell  and 
Thrush;^32^  (A),  Appleton  et  al;^38^  (B),  Byron;^3^  (C) ,  Cary.^33^ 
Full  and  dashed  curves:  see  text,  Sec.  3.2. 

Comparison  between  theoretical  predictions  of  the  recombination 
rate  coefficient  and  experimental  measurements  for  0^  as  a 
function  of  T.  •,  Campbell  and  Thrush;^37^  0,  Wray;^39^  (W) , 
Wray;^1^  (WM) ,  Watt  and  Myerson.  Dashed  curves:  see  text, 

Sec.  3.3. 

Comparison  between  theoretical  predictions  of  the  dissociation 
rate  coefficient  and  experimental  measurements  for  0^,  as  a 
function  of  T.  (W) ,  Wray;^^  (CV),  Camac  and  Vaughan; ^2^  (A), 
Anderson ;^3^  (WM),  Watt  and  Myerson. Dashed  curves:  see 
text.  Sec.  3.3. 

Comparison  between  theoretical  predictions  of  the  dissociation 

rate  coefficient  for  F^,  and  the  experimental  measurements 

(44) 

taken  from  Johnson  and  Britton.  Dashed  curve:  see  text, 


Sec.  3.4. 
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Figure  6  Comparison  between  theoretical  predictions  of  the  re* 

combination  rate  coefficient  and  experimental  measurements  for 

(45) 

C&2 ,  as  a  function  of  T.  §,  Clyne  and  Stedman;  (D) ,  Dlesen 

and  Felmlee;^6^  (C),  Van  Thiel  et  al,^^  Carabetta  and  Palmer, 

(49) 

and  Jacobs  and  Gledt.  Full  and  dashed  curves:  see  text,  Sec. 


3. A. 


Figure  7 


Figure  8 


Figure  9 


Figure  10 


Comparison  between  theoretical  predictions  of  the  recombination 
rate  coefficient  and  experimental  measurements  for  Br^,  as  a 
function  of  T.  t,  Ip  and  Burns;  (A),  Johnson  and 

Britton,  Warshay,^^  Britton,  and  Palmer  and  Homig,^8\ 
as  summarized  by  Ip  and  Burns;  (E) ,  emission  measurements  of 

(52) 

Boyd  et  al.  Full  and  dashed  curves:  see  text.  Sec.  3.4. 
Comparison  between  theoretical  predictions  of  the  recombination 
rate  coefficient  and  experimental  measurements  for  I2,  as  a 

(5g) 

function  of  T.  0,  Porter  and  Smith;  (BD) ,  Bunker  and 
Davidson; (BDGS) ,  Britton  et  al;^59^  (T) ,  Troe  and  Wagner. 

Full  and  dashed  curves:  see  text,  Sec.  3.4. 

Comparison  between  theoretical  predictions  of  the  dissociation 
rate  coefficient  and  experimental  measurements  for  HF,  as  a 
function  of  T.  0,  Jacobs  et  al;^8^  (B) ,  Blauer.^8^  Long 
dashed  and  short  dashed  curves:  see  text,  Sec.  3.5. 

Comparison  between  theoretical  predictions  of  the  dissociation 
rate  coefficient  and  experimental  measurements  for  HCi.,  as  a 
function  of  T.  Point  measurements:  Seery  et  al;^8^  (F), 

Fishburne ;  (J) ,  Jacobs  et  al.^8^  Long  dashed  and  short 

dashed  curves:  see  text,  Sec.  3.5. 


Figure  11  Comparison  between  theoretical  predictions  of  the  dissociation 


rate  coefficient  and  experimental  measurements  for  CO,  as  a 

(69) 

function  of  T.  Point  measurements:  Appleton  et  al; '  Shaded 
area:  Davies,  and  Presley.  Dashed  curves:  see  text, 

Sec .  3.6. 

Figure  12  Comparison  between  theoretical  predictions  of  the  recombination 
rate  coefficient  and  experimental  measurements  for  NO,  as  a 
function  of  T.  •,  Campbell  and  Thrush (W) ,  Wray  and  Teare.^^ 
Dashed  curves:  see  text,  Sec.  3.6. 
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